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Probability distribution for Q in open-universe inflation 

Alexander Vilenkin and Serge Winitzki 

(February 1, 2008) 
The problem of making predictions in eternally inflating universe that 

thermalizes by bubble nucleation is considered. A recently introduced reg- 

ularization procedure is applied to find the probability distribution for the 

ensemble of thermalized bubbles. The resulting probabilities are shown to 

be independent of the choice of the time parametrization. This formalism is 

applied to models of open "hybrid" inflation with O < 1. Depending on the 



parameters of the model, the probability distribution for is found to have a 

m 

peak either very close to = 1, or at an intermediate value of O in the range 



0.03 < n < 1. 



in 

O : I. INTRODUCTION 

OS 

A flat universe with fl = 1 is usually regarded as a firm prediction of inflationary models 
]]]]. However, observational evidence for a flat universe is far from being certain, and pro- 
gressively more precise estimates of the matter density of the Universe may yet show that 



the value of Q is less than 1. Not surprisingly, theorists found ways of modifying the models 
to make them compatible with Q < 1 |2|-|6|]. In the new class of models, called "open uni- 
verse inflation", the inflaton potential has a metastable minimum separated from the true 
vacuum by a potential barrier. The false vacuum decays through bubble nucleation, and 
the inflaton field rolls towards the true vacuum inside the bubbles, while inflation continues 
outside. Co-moving observers inside a bubble would, after thermalization of the inflaton, see 
themselves in an open homogeneous universe with Q < 1. In the usual inflationary models, 
flatness (Q ~ 1) results from a large amount of inflation needed for the observed homogene- 
ity of the universe. However, in models based on bubble nucleation, the homogeneity of the 
open universe inside a bubble is ensured by the symmetry of the bubble. So, the number 
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of e-foldings of inflation after nucleation (typically, of order 60) can be fine-tuned to give a 
specific value of Q between and 1. 

Although the nucleating bubbles expand at speeds approaching the speed of light, the 
false vacuum regions that separate them expand even faster. As a result, inflation never ends 
and bubble nucleation continues ad infinitum. If all nucleated bubbles are identical, then 
they will evolve to thermalized regions with the same value of Q. (To compare the values of 
Q in different regions, we can evaluate them at a fixed reference temperature, say, T = 2.7K.) 
However, there may be several types of bubbles giving rise to different values of Q, and in 
some models, like the model of hybrid inflation considered by Linde and Mezhlumian ||, Q 
can be a continuous variable. A natural problem in this kind of model would be to find the 
probability distribution for Q. This problem is the focus of the present paper. 

Our approach will be based on the assumption that we are "typical" among the civiliza- 
tions inhabiting the universe. Here, the "universe" is understood as the entire spacetime; our 
civilization is assumed to be typical among all civilizations, including those that no longer 
exist and those that will appear in the future. The assumption of being typical was called 
the "principle of mediocrity" in Ref. |7]]. It is a version of the "anthropic principle" which 
has been extensively discussed in the literature [|5|-|T2]]. In this approach, the probability for 
us to observe a certain value of Q is proportional to the total number of civilizations that 
will observe it. The total number of civilizations in a co-moving region can be expressed as 
the volume of that region V* at thermalization multiplied by the number z/ civ of civilizations 
that evolve per unit of thermalized volume [pT3| . The ratio of probabilities for a "typical" 



observer to find oneself in regions of type 1 and type 2 is then given by 



(1) 



T>(2) v (2) (2) • 

At this point the reader may be inclined to put this paper aside. What hope can we have 
to estimate the number of civilizations if we do not understand the conditions necessary for 
the evolution of life, let along consciousness? However, we believe that the situation is not as 
bad as it may seem. In models with a continuous spectrum of Q, like the model of Linde and 



Mezhlumian, the nucleated bubbles have identical particle physics. The difference in z/ C j v (Q) 
is then due only to the difference in the evolution of density fluctuations in bubbles with 
different values of Q. Roughly speaking, z/ civ (fl) is proportional to the density of galaxies 
formed in bubbles with the corresponding value of Q, and its calculation does not require 
any input from biochemistry. We shall attempt to estimate u civ (O) in Sec. |VT|, but until 
then our main goal is to develop a method for evaluating the volume ratios vi 1 ^ /vi 2 \ 

The volume V*^ 1 is the combined 3-volume of the hypersurfaces of constant temperature 
T = T* inside bubbles of type j (see Fig. [l|). These thermalization hypersurfaces are spaces 
of (approximately) constant negative curvature, and thus have infinite volume. In order to 
define the volume ratio in Eq. ([!]), this infinity has to be regularized. 

The most straightforward approach to regularization is to include in V* only the part of 
the volume that thermalized prior to some time, t = t c , and take the limit of the volume 
ratios V* /vi 2 "* as t c — > oo. However, the result of this procedure is highly sensitive to the 
choice of time variable t (that is, to the choice of the cutoff surface) flI^ , |T5|j . An alternative 
regularization prescription |TB[ is to cut off the volumes V* at the times te at which 



a small fraction e of the corresponding co-moving volumes is still in the inflating region. 
The probability ratios (|Ij) are then defined by taking the limit of V* /V*^ as e — > 0. The 
application of this procedure (which we shall call the e-prescription) to models of stochastic 
inflation was discussed in Refs. []16|,|17| , where it was shown that the resulting probabilities 



are essentially independent of time parametrization fT8|j . 

In models of stochastic inflation, the inflaton field undergoes quantum fluctuations on 
the horizon scale, and its evolution is described by a diffusion equation |19|] . The physics of 
such models is very different from that of bubble nucleation and expansion, and the methods 



of [|I6| , |T7|| are not directly applicable to open-universe inflation. The purpose of this paper 
is to extend the results of [16|;[T7| to this case and, as an application, to find the probability 
distribution for Q in hybrid inflation models. 

The outline of the paper is as follows. In Sec. [TJ we develop the geometric formalism 
necessary to describe the thermalization hypersurfaces within expanding bubbles. In Sec. 



I ll we calculate the regularized volume ratios in models with a discrete set of bubble types. 



Then in Sec. [TV] we verify that the volume ratios thus obtained are independent of time 



parametrization. In Sec. [V] we extend the analysis to the hybrid inflation model of || with a 
continuous family of bubbles. In Sec. [VT] we estimate the "human factor" v c \ w and calculate 
the probability distribution for Q. Conclusions follow in Sec. |V11| . Some calculations for 
Sections || and [nj are presented in Appendices. 

II. BUBBLE GEOMETRY 

The goal of this Section is to find the 3-volume of a thermalization hypersurface cut off 
at a given time t e . For simplicity, we shall use proper time for calculations; it will be shown 



in Sec. |V| that the resulting probabilities do not depend on the choice of the time variable. 

In models of open-universe inflation, the inflaton potential V ((p) has a local minimum, 
V — Vo, corresponding to a metastable false vacuum. In regions occupied by the false 
vacuum, the metric is approximately de Sitter, 

ds 2 = -dt 2 + exp (2H t) [dr 2 + r W] , (2) 

where dfl 2 = d6 2 + sin 2 9d(j) 2 is the usual spherical surface element, H is determined by the 



false vacuum energy, H = a/87tVo/3, and we use the Planck units, h = c = G = 1. 

At the moment of nucleation, a spherical bubble is formed, with the inflaton field in its 
interior on the other side of the potential barrier (with respect to the false vacuum). The 
bubble then expands and the inflaton field inside it evolves toward the true vacuum value, 
where it thermalizes. The interior of the nucleated bubble looks like an open Robertson- 
Walker (RW) universe in suitable coordinates (r, £), 

ds 2 = -dr 2 + a 2 (r) [d£ 2 + sinh 2 £dVl 2 ] . (3) 

The scale factor a (r) can be found from Einstein and scalar field equations, which in the 
slow roll approximation take the form 
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Equations (Hal)- 



and 



1 da 
a dr 



1 

' ^ 
dip 
d^ 



(4a) 
(4b) 



are valid provided that 
H' 



2nH 



< 1 



(5) 



Iff' I > H 2 . 



(6) 



Eq. (|5|) is the condition of slow roll, and Eq. (|) ensures that quantum fluctuations are 
small, so that the evolution of a and ip is essentially deterministic. The coordinates r and 
t can be chosen so that the center of the space-time symmetry of the bubble corresponds 
to t = r = 0. Then the surface r = is the future light cone of that center (see Fig. 0). 
We assume that the initial bubble size is small on the horizon scale Hq , so that for our 
purposes the boundary of the bubble can be approximated by this light cone. 

The relation between the coordinates (t,r) and (r, £) can be easily found if we assume, 
following || , that (i) the potential V (ip) has nearly the same value on the two sides of the 
barrier, and (ii) that the gravitational effect of the bubble wall is negligible. (A similar, 
although more cumbersome, calculation can be done for the more general case of non- 
negligible bubble wall gravity and different expansion rates H , Hi at the two sides of the 
wall; see below and Appendix B for details.) Then, at sufficiently small values of r, the 
geometry inside the bubble is close to that of de Sitter space with the expansion rate H . 
The solution of Eq. (ga|) with H {<p (r)) w H is 

1 



or 



sinh (H t) 



(7) 



This is accurate as long as 

\AH\r 



dH dip 



dip dr 
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(8) 



which gives 

r^lH'Mf 1 . (9) 

Here, ipo is the value of the field immediately after tunneling. At times r satisfying ([]) the 
coordinates (r, £) are related to (t, r) by the usual transformation between spatially flat and 
open de Sitter coordinates: 

t (r, £) = — In (cosh H t + sinh H t cosh £ ) , (10a) 
Ho 

r (r, = — sinh g r sinh e (1Qb) 

i/o cosh i^o'r + sinh f/oT cosh £ 

For r > jiif (</? )|~ , Eqs. fllCf ) no longer apply, but the coordinates (t,r) can be continued 
to the entire bubble interior as co-moving coordinates along the geodesies r = const. 

Thermalization occurs at a hypersurface of equal RW time, r = r*. The time r* of 
thermalization is found from the evolution equation fl4b"|) : 

where y?* is the value corresponding to the end of the slow roll regime near the true vacuum. 
We shall assume that i/o 7 "* 3> 1- 

The cutoff of the thermalization hypersurface at a time t = t e corresponds, in terms of 
the RW coordinates (r, £), to cutting off the surface r = r* at some value £ = £*, where £* is 
found from the requirement that the proper time t at (r*,^) be equal to t t . Therefore, we 
need to find the proper time along a geodesic r = const which starts in the false vacuum, 
continues into the bubble, and ends at (r*,£*). This task is facilitated by the observation 
that the time t along a co-moving geodesic in the de Sitter space after crossing the bubble 
boundary (for exp (H r) ^> 1) becomes almost identical to the RW time r inside the bubble: 

t (r, = r + A In cosh | + O (e- HoT ) . (12) 

In Appendix A it is shown that a co- moving geodesic r = ro, after crossing the bubble, 
rapidly approaches the RW co-moving geodesic line £ (r) = const. This, as well as Eq. 



(|12|), holds for times r within the range (^|), when deviations of the bubble interior from 
de Sitter space are small. But since the geodesies r = const and £ = const nearly coincide 
at r > Hq 1 , it is easily understood that Eq. (|1^) is valid throughout the bubble interior. 
Hence, the condition for the cutoff becomes: 

2 A, 

t* + — In cosh = t e . (13) 
-"0 2 

The solution of (|13|) can be written as 

= 2 cosh' 1 exp Fo(t ' 2 ~ r * ) . (14) 

The part of the thermalization hypersurface r = r* we are interested in is bounded by 
< £ < £*■ Its 3- volume, calculated using the metric (0), is 

K (t e ) = 4vr a 3 (tv) sinh 2 £d£ = ™ 3 (r*) (sinh 2£* - 2^) . (15) 

The calculations of the time cutoff ( |i~3| ) were performed for the case of unchanged ex- 
pansion rate Hi = H Q in the bubble interior immediately after nucleation. The analogous 
cutoff condition for H\ ^ Hq is derived in Appendix B. The thermalized volume V* (t e ) as a 
function of is still given by fllpp. 

III. REGULARIZED VOLUME RATIOS 

In this Section, we consider the situation where the bubbles come in several types (labeled 

by 1, 2, etc.). We assume that the nucleation rates for bubbles of each type are 71, 72, etc. 

A straightforward generalization to a continuous variety of bubbles will follow in Sec. 0. 

Our purpose is to find the thermalized volume ratios in bubbles of different types. For that, 

(1 2) 

we need to find the cutoff times te ' and evaluate the ratio of volumes of thermalization 

(1 2) 

hypersurfaces regularized by cutoffs at t e . To simplify our calculations, we shall first 
consider nucleation of bubbles of one type with nucleation rate 7, and subsequently generalize 
to multiple types. 
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For a bubble that nucleates at time t = 0, the regularized volume of thermalization 
hypersurface is given by Eq. ([15]) of the previous section. Now we have to account for 
bubbles nucleated at all times, starting for convenience at t = 0. Bubbles will nucleate 
in spacetime regions that are not already inside bubbles. (We disregard the possibility of 
tunneling from the true vacuum back to the false vacuum.) A point (to, r) will not be inside 
a bubble if no bubbles were formed in its past lightcone. The volume of the past lightcone 
of the point (to, r = 0) in de Sitter spacetime is 

f t0 47T 

V ic (to) = —rf c (t)exp(3H t)dt, (16) 

where r\ c (t) is the null geodesic ending at time to at r = 0, 

ri ^ = exp (-H t) - exp (-H t ) 
Ho 

This gives 

Vi c (to) = ^| - + O (exp (-H t )) . (18) 

Therefore, for sufficiently late times t 3> Hq 1 , the probability for a point (t,r) not to be 
inside a bubble is 

/ 4-7T7 \ 

^outside (t) = exp (-jVic) = exp y-^jptj , (19) 

where we assumed that the nucleation rate is small, 

l/H 4 « 1, (20) 

and accordingly disregarded the factor exp (22j/9Hq). 

The cutoff time t e is found from the condition that a fraction e of co-moving volume is 
still inflating at that time. Since inflation continues for some time inside the bubbles, the 
probability P; n f (t) of a point (t, r) to be in a still inflating region is not the same as the 
probability ([Lj]) of being outside bubbles. If we assume that inflation inside bubbles lasts for 
a period of proper time approximately equal to (the thermalization time given by (|TTD ), 
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then the points that are still inflating at time t are those which were outside bubbles at time 

t - r*: 

Pnf (t) « Poutside (t ~ n) . (21) 

Eq. ( PT| ) is not exact because the proper time t is different from the time r measured by 
the co-moving clocks inside the bubble; however, this difference is not large because the 
co-moving geodesies that define t quickly approach the RW geodesies inside the bubbles 
soon after they cross the boundaries. We will show in Appendix A that Eq. fl2l|) is accurate 
as long as the nucleation rate is small as assumed in (|20|). 
Hence, the cutoff condition becomes 

Pint it) = P outside (t e - r*) = e. (22) 

In (|22"D, we can use the asymptotic formula (|19|) for P ou tside (t) because we will be taking the 
limit of e — > oo for which H Q (t e — r#) ^> 1. 

Consider now a co-moving spatial volume equal to CHq 3 at t = 0, where C is a nor- 
malization constant corresponding to the initial number of horizon-size regions. The total 
volume of regions outside bubbles at a later time t 3> Hq 1 is given by 

Outside (*) = CH 3 exp (3H t) P ou tsidc (t) w CP - 3 exp (rf^ t) , (23) 
where d is the fractal dimension of the inflationary domain, 

^ 3 -§- < 24 > 

We will later use the fact that d ~ 3. 

The number of bubbles nucleated within the time interval (ti,ii + (iti) is 

diV(ti) = TKutside (*i) (25) 

and therefore the combined thermalized volume inside all bubbles (from t = until the 
cutoff time t f ) is 



V* = / K (t e - t x ) dN {h) = [ K (t e - ti) 7Kutside (ti) dti, (26) 
Jo Jo 

where V* (t) is given by (|l^). The integration in ( j2~6"|) is until t e — r* because bubbles nucleated 



after that time will not thermalize before t e . We can use the asymptotic formula ([23|) for 
Kutside (t) since the integral in (|26|) is exponentially dominated by bubbles nucleated at late 
times. 

Substituting (|T3]) , (|T5D and (p3| ) into (p6|), we obtain: 



V* = C—^a 3 (r*) / exp (cf# *i) (sinh 2£* (i e - * x ) - 2£* (t £ - h)) dh 

= C-^a 3 (r*) exp [d# (*e - r*)] / exp (-dtf t P (0) 2£ - 2£) -^£. (27) 



Here, £* (t) is the solution of fll3|) with t instead of £ e at the right hand side, and £ p (£) is the 
inverse function, 

ipK) = l-l„cosh| = l--A.n T A zf . (28) 

The time £ p (£) + is the proper time until thermalization along a co-moving geodesic 
that thermalizes at a given value of £; the formula (P8[ ) was derived for the simple case of 
unchanging expansion rate H . In Appendix B we find, for the case of Hi/H = h ^ 1, an 
expression for t p (£) similar to (f28|) : 

„ otp(0 ^_i±A ln _i±l_ (29) 

This coincides with (]2q) for h = 1. 

The integration in (|27D is performed up to £ max = £* (t e — r*). Since d ~ 3, and (0 ~ 
£ for £ ^> 1, the integrand of (p7|) decays exponentially at large £, so the precise value of 
£ max is unimportant, and we can take the limit £ max — > oo. The resulting integral with i p (£) 
given by ( |29|) depends only on h = Hi/H and can be expanded in (3 — d) as 

r Si ^W^ o^ = f(h)+0(3-d), (30) 
where the function / (h) can be approximated |2(J within an error of 2% by 
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(3D 

Keeping only the leading term of the expansion in (3 — d), Eq. fl27|) for the thermalized 
volume becomes 

V* = Cf (h) ^ exp (dH (t e - rO) al (32) 
-"o 

where a* = a (r*). 

The cutoff time £ e is found from (p2|) , 

exp [- (3 - d) # (^ - r0] = e, (33) 

and we obtain, after substituting in (^) and simplifying, 

V* = Cf(h)^e-^al (34) 

The expression fl34|) for the thermalized volume holds if there is only one type of bubbles. 
In the case of several bubble types, the argument above is modified in the following points: 
(i) the nucleation rates 7^ , the thermalization times r* and the volume expansion factors 
a* are specific for the j-th type of bubbles; (ii) the fractal structure of the region outside 
bubbles is affected by nucleation of bubbles of all types; the corresponding fractal dimension 
d is 

-'=»-^E^-3-^ (35, 

(iii) the cutoff condition (^) is modified for bubbles of type j to 

^outside (t® - t&) = exp (- (3 - d) H (t?) - r^))) = e. (36) 

The motivation for (|36|) is as follows. The cutoff procedure for bubbles of type j sets the 
cutoff time te at which a fraction e of all co-moving volume that will eventually thermalize 
in bubbles of type j, is still not thermalized. Since bubbles nucleate at time-independent 
rates 7^ per spacetime volume, the probability for a given observer outside any bubbles to 
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thermalize in a bubble of type j is at all times proportional to 7^). Therefore, at any time 
t, a fraction 7^') /y of the co- moving volume that is outside bubbles at time t, and the same 
fraction 7W/7 of the total co-moving volume, will eventually thermalize in bubbles of type 
j. According to Eq. (fL9D , a fraction exp (— yV\ c (£)) of all co-moving volume is still outside 
bubbles at a time t; then also a fraction exp (— yV\ c (t)) of the co-moving volume that is to 
thermalize in bubbles of type j, is outside bubbles at time t, and this holds independent of 
j. Hence the cutoff condition fl3"3|) is only modified for a given type j in its dependence on 
d and t*, as written in 



The regularized thermalized volume V* corresponding to bubbles of type j becomes 

(?) 

V« = C7/(/iW)^e-^[a?)] 3 . (37) 
The ratio of volumes in bubbles of types, e.g., 1 and 2 is 



VP_ yW 

V (2) « 7 (2) 



al 2) 



3 



/(,!<1)) (38) 
Since the ratio is independent of e, the ratios of thermalized volumes in bubbles of different 



types are directly given by Eq. (B9) |21|. 



IV. ARBITRARY TIME VARIABLES 

We consider now a different choice of time variable i related to the proper time t, along 
a geodesic r = r , by: 

dt = T(H(t,r ))dt, (39) 

where T (H) is an arbitrary (positive) function. Such a relation will, for instance, describe 
the proper time (T = 1) and the "scale factor" time (T (H) = H). We can always normalize 
i so that T (H Q ) = 1. Then, the new time variable i will be identical to t in de Sitter regions 
where H = Hq. However, inside bubbles the time variable will be significantly changed. In 
this Section, we will modify the calculations of the preceding sections to accommodate the 
new time variable and show that the result (|38|) is independent of the choice of T (H). 
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As in Sec. ||, we calculate the time along a co-moving de Sitter geodesic by matching 
it with a Robertson- Walker geodesic at a time r . The thermalization time (|TT| ) is then 
modified to 

r* = t + I" *T(H(r))dT. (40) 

J TO 

The calculations of the co-moving and physical volumes outside of bubbles (|T^) — (^) and 



of the number of nucleated bubbles (25) concern only the de Sitter region, therefore for 
the new time variable the same expressions hold, and the fractal dimension d is unchanged. 
Equation (|22D for the cutoff t e is modified to 



outside 



(t e - f*) = e. (41) 



In the calculation of the thermalized volume (|T5|), the integration is performed on the ther- 
malization surface that does not depend on time parametrization, so the result (|15"|) holds. 
The spatial cutoff becomes 

£ = 2 cosh- 1 exp g°(jLr?») . (42) 

The regularized thermalized volume is found analogously to (|27|), except that the integration 
is done over the time of bubble nucleation i\ in the new time parametrization. The calcu- 
lations are identical, except for the changed values of f*, and the results (PT|), fl3"B|) depend 
on f* only through invariant factors a* (f*) given by 



Tt 



a* (t*) = exp / H (r) dr = exp 



o 



(43) 



where y9 and tp* are appropriate initial and final field values. We conclude that the regu- 
larized probability ratios (^) are independent of time parametrization. 

V. THE LINDE-MEZHLUMIAN MODEL 

Linde and Mezhlumian considered a model of hybrid inflation in which homogeneous 
open universes with different values of Q < 1 are created via bubble nucleation. In that 
model, two scalar fields a and <fi evolve in an effective potential of the form 
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V{a, ( f>) = V (a) + (T 2 V 1 ( ( f>), (44) 

where the potential Vq (a) has two minima corresponding to the false and true vacua, re- 
spectively (Fig. |IJ), and V\ ((f)) is some potential with a slow-roll region suitable for "chaotic" 
or "new" inflation. While the field a stays in the false vacuum (a = 0), the potential for 
is flat, and quantum fluctuations smooth out the distribution of to almost uniform (up to 
corrections due to tunneling, see below). To make this distribution normalizable, we shall 
assume that is a cyclic variable and identify = <pc with = 0. The field a has a small 
probability to tunnel to the true vacuum through the formation of bubbles which will have 
a continuous spectrum of values of 0. Inside a bubble, the potential becomes 0-dependent 
and the field starts evolving from its initial value 0o until thermalization in the global 
minimum of V (a, 0). Depending on the initial value 0o, the bubbles will undergo different 
amounts of inflation and, therefore, will have different values of Q. We shall apply the re- 
sults of Sec. H| to calculate the probability distribution for Q in this ensemble of bubbles, 
for a particular family ( |44"D of potentials V (a, 0). For sufficiently large values of V\ (0) the 
tunneling is absent, since the a 2 term raises the true vacuum energy above that of the false 
vacuum. We can choose the potential V\ (0) so that tunneling is allowed only for values of 
satisfying ([]), and thus quantum fluctuations of will not be dynamically important inside 
the bubbles. We shall also assume that the value of does not change appreciably during 
tunneling. 

The type of bubble is now characterized by a continuous parameter O , the value of 
at tunneling. To apply the result of Sec. |T|, we need to supply a measure in the parameter 
space, i.e. a weight for the bubbles with in the interval (0 O , 0o + d(p ). The situation differs 



from Sec. [Ill] also in that the nucleation of bubbles of different types occurs in different 
regions of space. To account for this, we describe the inflating regions of false vacuum by a 
stationary solution of the diffusion equation for the volume P (0 O , t) d(p of regions occupied 



by the field in the interval (0o, 0o + <^0o) at time t ||19|| . The diffusion equation is modified 
to include a "decay" term for bubble nucleation: 
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47T7 (0) 

3H 3 



(45) 



Here, if (0) is the expansion rate in the false vacuum, and 7 (0) is the ^-dependent tunneling 
rate. For the potential fl4*4]), which is our concern here, H (0) = H = const. The stationary 
solution of (|45| ) can be written as 

P{<f>,t) = P {<f>)e W {dH t), (46) 

where P$ {(f)) is the highest eigenvalue solution of the stationary diffusion equation 

w t f 3 -^m Po=dPo (47 ) 



8tt 2 <90 2 V 3# 4 7 

with periodic boundary conditions, and d is the corresponding eigenvalue. According to Eqs. 

— (il)> the resulting thermalized volume in bubbles of a given type is proportional to 
the volume of the regions of false vacuum in which bubbles of that type can nucleate. The 



latter volume is proportional to Pq (0o) <^0o- Therefore, the probabilities of Sec. |TTT| should 
be weighted with P (0 O ) <^0o- 

By integrating Eq. (|4~7|) over 0, we obtain an expression for d: 

i = 3 _*^SfMM* (48) 

Since the tunneling rate 7 is small, we can approximate the solution of ( [4TD by a constant 
function, and then the eigenvalue d is given by the formula similar to (j35|): 



4tt 1 



<t>c 



^^3-— i — / 7 (</>)#■ (49) 
3H^ (j) C Jo 



According to (|38|), the probability distribution depends on O through the nucleation rate 
7 (0o), the expansion factor a* (0o), an d the factor / [H (0 O ) /Ho] = f (0o) which describes 
the effect of a different expansion rate H\ = H (0o) i n the bubble interior after nucleation. 
The nucleation rate per unit spacetime volume is estimated || using the Euclidean 0(4)- 
symmetric instanton solution o (r) for the field a coupled to gravity: 
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7 (0o) = A(0 Q )exp(-£ E (0 o )), (50) 

where Se (0o) is the instanton action and A (0 O ) is the prefactor which we assume to be a 
slowly-varying function of 0o- The regularized probability of being in bubbles that tunneled 
with = 0o is then expressed, with a suitable normalization constant N, as 

c£P(0o) =iVz/ civ (0 o )rfP(0 o ), (51) 

where we have separated the distribution dV (0 O ) due to the thermalized volume, 

dV (0 O ) = A (0 O ) exp (-S E (0o)) 4 (0 O ) / (0o) (0o) d0 o , (52) 

from the "human factor" z/ C i v (0o) introduced in Eq. (|l|). We will now concentrate on the 
above distribution, whereas the effect of the factor z/ civ (0 O ) will be discussed in the next 
Section. We shall be interested in the leading (exponential) dependence on O i n Eq. (|52|) 
and shall therefore approximate the factors A (0o), f (0o) and Pq (0o) by a constant. 

The expansion factor at thermalization a* (0o) for a bubble formed at the value = 0o 
is determined by ([43|), 

'+* V(a ,<f>) J 



a* (0o) = exp 



-87T 



'<f>o K^l^O,' 

where the value 0* corresponds to the end of slow roll and is defined by 

V 



(53) 



4ttV 



1. (54) 



Eqs. (il) — (^1) give the probability distribution for the value 0o at which tunneling of 
the field a occurs. To obtain a probability distribution for Q, we need to find the present 
value of Q as a function of 0o- As outlined in [[|, we can relate Q to the expansion factor 



-D \ - I J-thX J-eq 



a* (0o) given by ([53]): 

^(0o)=fl + ^- T ) , 5=1^1 =^_, (55) 

V «*(0o)/ V 7 "// T CMB 

where is the thermalization temperature, Tcmb is the cosmic microwave background 
temperature at present, and T e(? is the temperature at equal matter and radiation density. 
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Depending on T t h, the value of B is ~ 10 25 — 10 50 . A higher value of V\ (0o) corresponds to 
longer inflation and a larger expansion factor a*, and therefore to a value of Q closer to 1. 

To calculate dV (Q) /dQ, we choose a potential that in the range of </> where tunneling is 
allowed is given by 

V(a^) = V (a) + ^ 2 a 2 , (56) 

where Vo (cr) still has the shape shown in Fig. [| A similar potential was also considered in 
f|. Note that the slow roll condition @ requires ^> 1. To facilitate the calculation of the 
instanton action Se (4>o), we shall choose Vq (a) to be quartic in a: 

V Q (cr) = Act 4 - b x a 3 + b 2 cr 2 + const. (57) 

The constant is chosen so that the true vacuum energy is zero, giving a vanishing cosmolog- 
ical constant. Since we assumed that the bubble size is small on the horizon scale, we can 
disregard the effect of gravity and treat the instanton as in flat space. The calculation of the 
instanton action in flat space for general quartic potentials of the form Q57D was performed 
semi-analytically in [E^], and we shall use the result obtained there, 



7T 2 a x 5 + a 2 S 2 + a 3 S 3 
M 0o) = 3A ^3 , (58) 

where a,\ = 13.832, a 2 = —10.819, a 3 = 2.0765 and the dimensionless parameter 5 (0 O ) is 
defined, in terms of the parameters of the potential (|57j), by 



5(0 O )=8A - \ 2V0 . (59) 



The allowed range of 5 is from its minimum value <5 m ; n = 8A62/&1 to 2, where 5 = 2 corre- 
sponds to the maximum value of <fi at which tunneling can still occur:, 

2 / h 2 \ h 2 
±2 z / °i 1. \ u i 

rmax 



- 9 {tx- b *)=Tt g {2 - s ^- < 60 » 

The thin wall approximation, valid when the minima of the potential fl57|) are almost 
degenerate, corresponds to <5 m i n ~ 2, and then a±5 + a 2 S 2 + a 3 5 3 ~ 1 for 5 min < 5 < 2. A 
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generic choice of parameters A, b\ and b 2 , such as b\ ~ Ao"o, b 2 ~ Acrg, will give <5 min ~ 1. 
Then, the expression a±6 + a 2 8 2 + a 3 5 3 in fl5"8|) is also of order 1 for the allowed range of S. 
Accordingly, we will disregard this expression below. 

Using the potential (0) — (0), we can calculate a* (0o): 



a* (0o) = exp 



S7T 



exp 



8tt 



00 | 2 a 2 + Vq ((Tq) 

0. £ W + ^oW)' 



(61) 



Since the potential (|56"D includes interaction between a and 0, the value of the true vacuum 
a will be slightly 0o-dependent, o"o = cr (0o), and as the field slowly evolves toward 
= 0*, the field a will follow the shifting position of the minimum o"o (0o)- Without the 
dependence of o"o on 0, the expansion factor would be 



a* (0o ) = exp 



57T 



exp (27T02 - 2tT0 2 ) . 



(62) 



The exact expression for a* (0o) contains a correction to (|62] 



a* (0o ) = exp 



\ \ rma> 



(63) 



where the function F behaves as F (x) ~ x 2 at small a; and F (1) < 1 (the explicit form of 
F is unimportant). 

Eq. (|55| ) for (0 O ) becomes 



CI 



[l + J Bexp(-47r(l + F)(02-02))] 1 



(64) 



Assuming that In B ^> 1, one can see that Q (0o) changes very quickly from to 1 in a narrow 
region of relative width A0/0 ~ (In I?) 1 around = 0i, where 0i = a/ (l/47r) hi I? + 2 . 
For a typical value of In B ~ 100, one obtains 0i ~ 3. Note that for 0o < 1 the slow roll 
approximation is not valid and Eqs. (|61~D — (|64"D are not applicable; we shall only consider 
the distribution ( |52"D for O > 0*, where 0* ~ 1. Correspondingly, the range of Q is from 
^ (0*) ~ B~ l « to fi max = ^ (0max)- The maximum value of Q is 



fimax = [1 + B exp (-47T (1 + F (1)) 0^)] 1 . 



(65) 
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Generically, max 1 and f2 max is very close to 1. 

Combining (0), @), and (|63|) , we obtain the leading exponential dependence of the 
distribution ( |5T| ) on </> : 

oc exp (-^ (2 - 5 (0 o )r 3 + 6vr (1 + F) - 2 )) . (66) 



VI. PROBABILITY DISTRIBUTION FOR Q 

To obtain the probability distribution for Q = Q (<po), we need to transform dV (0o) to 
the new variable Q via 

Expressed as a function of 0o, this distribution is 

Since most of the range of Q (except a narrow region around Q — 1) corresponds to 0o <C 
0max, we can expand (|68|) in ^o/^max and obtain an approximate power-law dependence 

oc fi 1 / 2 " 3 ^ (1 - fi) 3 ^ 5/2 , (69) 

where we have defined the dimensionless parameter [i by 

fx = - (70) 

12A0Lx (2 - 5 min ) 3 

For f2 very close to 1, the right hand side of ( j68|) is dominated by the first term in the 
exponential, which makes it rapidly drop to 0. Depending on the value of /i, there are three 
distinct behaviors of dV/dQ (Fig. [|). In the first case, /i < 1/6, the function monotonically 
grows with Q until it peaks at Q — f2 pca k ~ fi m ax ~ 1 and very rapidly falls off to 
for Q > fipeak- The second case occurs for fi > 5/6; the distribution (|6T| ) monotonically 
decreases with Q and its maximum is at the lower boundary Q = 0. Lastly, in the third 
case, with 1/6 < fi < 5/6, the distribution ( |68D decreases from a local maximum at Q = 
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and then increases to another local maximum at fi = fi pC ak ~ 1 (Fig. 0). To determine 
which maximum dominates the probability distribution, we consider its approximate form 
(|B5D. If fi < 1/2, the second exponent in Eq. fl69|) is smaller than the first one, giving a 



stronger peak at fi ~ 1. For p > 1/2, the peak at = is stronger [f23 . 

Now we consider the influence of the factor z/ civ (0 O ) on the distribution In our 

model, the low-energy physics is identical in all bubbles, and therefore u C i v is simply pro- 
portional to the number of potentially inhabitable stellar systems. The structure formation 
process in different bubbles is also essentially the same, apart from the difference in fi. The 
main effect of fi is to terminate the growth of density fluctuations at redshift 1 + zq ~ fi -1 



24J. Assuming that the dominant matter compontent is "cold", density fluctuations begin 
to grow at redshift of matter and radiation equality, 1 4- z C q ~ 2 x 10 4 fi/i 2 . With h = 0.7, 
the overall growth factor is 

/ (fi) = \^ - 10 4 fi 2 , (71) 

where we have assumed that z eq > zq, that is, fi > 10~ 2 . Otherwise, there is no growth, 
and thus / (fi) ~ 1 for fi < 10~ 2 . 

If density fluctuations are generated by inflation, then their initial amplitude on each 
scale has a Gaussian distribution. Its rms value at horizon crossing, (8p/p) vms = 5, is 
determined by the shape of the potential V (ip, a) and is approximately scale-independent 
on astrophysically relevant scales. In bubbles with values of fi such that / (fi) 5 > 1, most 
of the matter is captured into bound objects, and z/ c ; v is essentially independent of fi. On 
the other hand, if f(Q)S <C 1, then almost no structure is formed. In this case, bound 
objects are formed only in the rare regions where Sp/p exceeds the rms value 5 by a factor 
> [/ (fi) S\ 1. Hence, we expect that in the range 10~ 2 < fi fi, 

z/ civ (fi) oc exp (— Kfi -4 ) , (72) 

where k ~ lO" 8 ^ 2 and fi ~ 10' 2 5^ 2 is the solution of / (fi) 6 ~ 1. For fi < 10~ 2 , / (fi) ~ 1 
and we expect z/ civ (fi) oc exp (— k'5~ 2 ) with k' ~ 1. The function z/ civ (fi) is sketched in Fig. 
^| for the full range of fi [ 25 1 . 
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The effect of the factor (O) on the probability distribution ( |6^ ) can now be easily 
understood. If dV/dTL has a single peak near Q — 1, as in Fig. |]a, then the peak position 
remains essentially unchanged, and the distribution function is suppressed only for Q < Q, 
where it was already very small. The most interesting modification occurs when there is 
a (local) peak at Q = 0, as in Figs. [|b-c. This peak is then shifted to a larger value, 
fipeak ~ y/2(3n- l/2) _1/4 n ~ a (For 5 ~ lO" 3 , 5 ~ 1(T 2 , and 5 ~ 10 -1 , we obtain 
Q ~ 0.3, Q ~ 0.1, and Q ~ 0.03, respectively.) In the case of \l > 5/6, this is the only 
maximum of dV/dfl. The behavior of the full probability distribution fl5lD , 

dV (O) = Nu civ (Q) dV (Q) , (73) 

is sketched in Fig. |^ for all three cases. 

The idea that anthropic considerations make a low value of Q very unlikely has been 
previously discussed by a number of authors 0,0; however, to our knowledge, no attempt 
has been made to make this argument quantitative. A similar approach to the cosmological 
constant has been developed in Ref. [p7|,EM[7],E3 . 



VII. CONCLUSIONS 

In this paper, we have considered scenarios of open-universe inflation, where the 
metastable false vacuum decays by quantum tunneling and forms bubbles of different types. 
The interior of a bubble is observed as an open universe, in which inflation continues until 
thermalization. Our goal was to find a probability distribution for thermalization in dif- 
ferent bubble types, following the approach of [|7|,p!6|, according to which the probability is 
proportional to the number of civilizations that will evolve in bubbles of each type. The 
problem of calculating the probability splits into a calculation of the ratio of physical vol- 
umes thermalized in different types of bubbles and of the number of civilizations u civ that 
evolve per unit thermalized volume. 

In Sections we developed a method for calculating the volume ratios. We first 

considered the case of a discrete set of bubble types and then extended the analysis to a 
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continuous spectrum of bubbles. As an example of the latter, we focussed on the Linde- 
Mezhlumian model of hybrid inflation H which gives rise to an ensemble of open universes 
with Q < 1. Since all nucleating bubbles in this model have identical particle physics, we 
were able also to estimate the "human factor" z/ C j v (f2). We found that, depending on the 
dimensionless parameter /i defined in ([71]), the probability distribution dV/dQ is peaked 
either at Q = 1 (for fi < 1/6) or at an intermediate value Q = Q in the range 0.03 < Q < 1 
(for /i > 5/6). For 1/6 < < 5/6, the distribution has local maxima at both Q = Q and 
ft = 1; the relative magnitude of the two peaks depends on /i and on the amplitude 5 of 
density fluctuations. 

In this paper, we have only considered models in which the false-vacuum regions inflate 
at a constant rate H. In a more general situation, the field <fi would evolve before as well 
as after tunneling, leading to a slowly changing H. The analysis of such models would be 
substantially more complicated, while the results are likely to be similar to those in simpler 
models with a constant H. 

APPENDIX A: PROPER TIME IN THE BUBBLE INTERIOR 

Here it will be shown that a co-moving geodesic continued from the de Sitter region to 
the bubble interior exponentially approaches a Robertson- Walker (RW) stationary geodesic 
inside the bubble. We shall calculate the proper time along such a geodesic and show that 
the approximate formula (|2lD is accurate within our assumptions. 

As we noted in Sec. |I|, there is a region inside the bubble in which the spacetime is 
approximately de Sitter, and the RW coordinates (r, £) in that region are related to de 
Sitter ones by (pi]). The range of r in that region is r 1/ H 1 ((fo), as follows from @. We 
can use the coordinate change (|T0"| ) to continue a co-moving geodesic r = r from the false 
vacuum region to the bubble interior (provided that the geodesic intersects the bubble, i.e. 
that H r < 1). The resulting trajectory £ (r) is 
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H r cosh H t + J sinh 2 H q t + (H Q r ) 2 

^ T ) =ln n n"wuf • ( A1 ) 

(1 - iJ r ) smh H r 

At large values of r such that expH^T ^> 1, the trajectory ( |A1|) becomes 

i.e. it is exponentially close to the co-moving geodesic line £ = const in the RW region. 
We see from (^), (H) that there is a range of r such that 

1<<7/ » r<< TO (A3) 

and in this range the co-moving world-lines, r = r , continued from the region outside of the 
bubble into the interior, become very close to the RW co-moving world-lines, £ = £o, while 
the spacetime is still sufficiently close to de Sitter. At times r satisfying ( |A3| ) , the proper 
time along r = r becomes exponentially close to r, as shown by (p~|). 

Now we will consider Eq. (^Tj) which was based on the assumption that the time interval 
between crossing the bubble boundary and thermalization is equal to r* for all geodesies. 
This assumption is not exactly true, because during the time period when the time variables 
t and r differ significantly, their difference depends on the spatial coordinate tq, which 
varies among different geodesies r = r . As a result, the proper time interval along a 
geodesic between entering the bubble and the point (r, £ ) differs from r by an r -dependent 
correction At. Assume for simplicity that the bubble is centered at r = 0. A co-moving 
geodesic r = r entered the bubble at time t given by 

*o = -(l/#o)ln(l-#oro). (A4) 

The correction At is then 

At (r ) = * (r, £ (r )) - to - r = J- In cosh ^ - t = —j- In (1 + H r ) . (A5) 

The function At (r ) does not depend on r and its maximum value is —Hq 1 In 2 (for H r w 
1)- 
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Now we can show that the correction ( |A5|) does not significantly influence Eq. ([H]). A 
change in the thermalization time r* by the correction At ~ Hq 1 in ( [21] ) would change Eq. 
(pl|) by the factor exp (—Attj/ (3Hq)), which is very close to 1 because, as we assumed in 
(p0|), 7/ (3Hq) <C 1. Therefore, Eq. ( pT| ) is accurate within our assumptions. 

APPENDIX B: THE CASE OF DIFFERENT EXPANSION RATES Hi £ H 

Here we present the calculations of the proper time until thermalization in the general 
case when the gravitational effect of the bubble wall is not assumed to be small and the 
expansion rate Hi inside the bubble significantly different from Hq (presumably, Hi < Ho) . 
We shall assume, however, that the size of the nucleated bubbles is small on the horizon 
scale Hq . 

The de Sitter spacetime is represented by the hyperboloid 

C 2 + w 2 - = H^ (Bl) 

embedded in a 5-dimensional space ((j, (2, C3? w i v ) with Minkowskian signature. For simplic- 
ity, we treat the bubble interior also as a de Sitter spacetime region with constant expansion 
rate Hi. Then the bubble interior will be a piece of the hyperboloid 

C 2 + (w - Aw) 2 -v 2 = Hf (B2) 



cut out by intersection with ( pl|) . The displacement Aw is related to the bubble wall 
tension or, alternatively, to the initial bubble size ^Of. The bubble wall is at w = wq = 
(1/2) [Aw - (tff 2 - H 2 ) /Aw]. 

The flat RW coordinates (t, x) in the outer region are introduced by 

H t = lnF (w + v) , (B3a) 

# x = (B3b) 

w + v 

This gives the trajectory of the bubble wall in these coordinates, 
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H r (t) = ^(1- e- H ^f + 2e~ H ot (l _ H w ). (B4) 

The assumption of small initial bubble size corresponds to H w ~ 1, which means that we 
can approximate the bubble wall by the lightcone H ri c (t) — 1 — e~ Hot . This considerably 
simplifies the algebra. 

Our goal is to find the proper time until thermalization along a co-moving geodesic that 
starts as r = const in the outer region and crosses the bubble wall. We introduce the flat 
RW coordinates (t^Xx) also in the interior region: 

if i*i = In ff i (w-Aw + v), (B5a) 

ffixi = (B5b) 

w — Aw + v 

The two coordinate systems are matched at the bubble wall, and the metric is continuous 
across the wall. This allows us to continue the geodesic r = r through the bubble wall 
by requiring that the component of its 4-velocity parallel to the wall be continuous. We 
denote by (3 the r component of the initial 4-velocity, (3 — found from this condition. 
The (generally non-zero) velocity ^ with which the geodesic emerges in the interior is 
determined by f3. A general radial geodesic in the interior de Sitter region is described by 



H in (ii) = r 10 + y/p-* + e^'io - y/P~ 2 + e H ^, (B6) 

where (tio,r 10 ) is the initial point at the bubble wall in the coordinates (ii,Xi) and P is a 
constant of motion related to the initial velocity by 

exp 



*10 



1 - exp {2H l t w ) 



t 2 

dti 



/3exp(^ 10 ). (B7) 



tio / 

The proper time <5t along this geodesic from the bubble wall crossing until time t\ is found 
to be 



sinh" 1 -p- - sinh" 1 — — J 



(B8) 
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The geodesic (p6|) asymptotes to the line r = r a at large times, where r a is given by 

1 ^exp(-gitio) 
ra " ri ° + ^ l + v /l + ^ • (B9) 

As in Sec. 0, we introduce the open RW coordinates (r, £) in the interior and match the 
geodesic ( |H6| ) with a line £ = at time t\ (r , £*) given by ( |10a| ). This enables us to find the 
total proper time until thermalization as the sum of the time to until bubble wall crossing, 
the time 5t {t\ (r , £*)) from the wall crossing to matching with £ = £ , and the time r* — r 
until thermalization: 

ttotai (Co) =t + 6t (tx (t , Co)) +n- T . (BIO) 

The trajectory flB6|) is completely specified by its asymptotic value of £, and we can express 
the parameters f3, P, r 1; t w and t through £ . After some algebra, we arrive at the following 
expression for the time ( piO|) : 

If l + l + h \ /N 

*„„({„)= r. + -({„-— In TT7 - r ^j. (Bll) 

For h = 1, this reduces to the left hand side of (|T3|), as expected. 

In the calculation of the thermalized volume in Sec. |ITJ, we will use the function t p = 
^totai (£o) — r *j which has the meaning of the correction to the thermalization time: 



Again, for h = 1 this expression coincides with (|28|). 
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APPENDIX C: FIGURES 




FIG. 1. A conformal diagram of bubbles nucleating in an inflating background. The shaded 
regions of spacetime inside the bubbles are thermalized. The thermalization surfaces are the 
boundaries of these regions. They have an infinite 3- volume which can be regularized by introducing 
a cutoff hypersurface t = t c and keeping only the part of the volume below this hypersurface. 
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FIG. 2. Geometry of the bubble interior. 



V (o) 




FIG. 3. The shape of the potential V (a) in Eq. (44) 



31 



(a) (b) (c) 




FIG. 4. Probability distribution dV(fl)/dQ with </>max/^* = 100, shown logarithmically up to a 
normalization, (a): /U = 0.01. The peak at £1 ~ 1 is extremely sharp; the ratio of the values at 
Q = .99 and at = is ~ exp 20, while the peak value differs from that at O = by a factor of 
~ exp (25000). (b): jj, = 0.5, there are two local maxima near Q = and £1 = 1. (c): /U = 2. The 
function monotonically decreases. The maximum value at Q = differs from a typical intermediate 
value (O ~ 1/2) by a factor of ~ exp (60). 
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FIG. 5. Dependence of v c i v on f2 for 5 = 10 2 . The ratio of the maximum and the minimum 
values is ~ exp 5~ 2 . The origin on the vertical axis is arbitrary. 
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FIG. 6. Probability distribution dV/dVi, shown logarithmically up to a normalization, in the 
three cases corresponding to Fig. ||a-c. The parameter values are the same as those in Figs. ||, |5| 
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